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' The Bean- Livingston barrier at the interface of type-II superconductor/soft-magnet heterostruc- 

, tures is studied on the basis of the classical London approach. This shows a characteristic dependence 

' on the geometry of the particular structure and its interface as well as on the relative permeability 

of the involved magnetic constituent. The modification of the barrier by the presence of the magnet 
, can be significant, as demonstrated for a cylindrical superconducting filament covered with a coaxial 

^ ■ magnetic sheath. Using typical values of the relative permeability, the critical field of first pene- 

' tration of magnetic flux is predicted to be strongly enhanced, whereas the variation of the average 

critical current density with the external field is strongly depressed, in accord with the observations 
of recent experiments. 



PACS numbers: 74.25.0p; 74.25.Sv; 84.71. Mn; 41.20.Gz 

o 

O : I. INTRODUCTION 

U 

S^. Heterostructures on the macro or nano scales involving type-II superconductor and ferromagnet elements 
^ ' show great potential for improving superconductor properties such as critical currents and critical fields, 
and therefore have been extensively studied both experimentally and theoretically during the past few 
a , year,gii^i3,4,5,6,7,8,9,io,ii,i2,i3,i4^i5,i6,i7,i8,i9,20,2i,22,23,24,25,26,27^ jf hard magnets are used, the interaction of the mag- 

^ ' netic vortices of the superconductor with the magnetic moments of the ferromagnet may lead to an enhancement 
of the pinning of the vortices^i^i^ or to an increase of the critical field a^i^^ . Soft magnets, on the other hand, aid 
to amend superconductor performance by shielding the transport current self-induced magnetic field as well as the 
externally imposed magnetic fieldiLiSiSSiSSiSS. Superconductors encompassed with such materials exhibit enlarged 
Q • critical currents through wide ranges of the strength of an applied field, when in the critical stateSSiSi^iiSSiSSr, and 
even overcritical currents, when in the Meissner stateSi. The latter state persists until the shielding and/or transport 
currents, which push magnetic vortices into the superconductor bulk, overcome the Bean-Livingston barrier against 
entry of magnetic flux"^''; an impediment created by the (positive) Gibbs free energy of the vortices themselves. This 
I suggests the surmise that, due to the induced magnetization, the presence of a soft magnet may alter the character- 
istics of nucleation of a vortex at the superconductor/magnet interface as compared to nucleation at the surface of a 
T— I superconductor facing vacuum, the phenomenon analysed hitherto22*2ii224^i24i2^. 

' Any theoretical study of the Bean-Livingston barrier at the interface of a type-II superconductor/soft-magnet 
[ heterostructure, discerned by the observable critical field of first penetration of magnetic flux and the observable 
^0 ■ average critical current density of loss-free transport of electric charge interlinked with it, must resolve two cardinal 
points: 

■ (a) the dependence of these observables on the geometry of the particular structure and its interface; 
(b) the effect of the relative permeability of the involved magnetic constituent. 

Here, we exemplify both traits for an infinite flat and, respectively, finite curved geometry of a type-II superconductor 
' next to a soft-magnet environment, adopting the classical London approach. 
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■ II. THEORY 
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\ The magnetic induction B in the superconductor region around a vortex obeys the London equation2& 

B + A'V X (V X B) = Q, (1) 
with the London penetration depth A and the source Q at position r given by 



Q(r) = $oy"ds,5(s-r), (2) 



where $0 denotes the quantum of magnetic flux and 6 means the Dirac delta function, the integration extending along 
the vortex core. The magnetic field H in the magnet region and the magnetic induction B in the entire space satisfy 



2 



the Maxwell equations 

V X H = and V • H = 0. (3) 

To simplify later analysis, we postulate the relationship B = /i/igH in the region confined to the magnet itself, 
assuming a (field-independent) relative permeability of the magnet, fi, apart from the permeability of free space, fi^. 
Furthermore, to avoid complications arising from the proximity effect, we invoke the presence of an insulating layer 
at the superconductor/magnet interface, of thickness much smaller than the London penetration depth (as observed, 
e.g., in MgB2/Fe composites2&) regarding, for mathematical convenience, this layer as infinitesimally thin. Boundary 
conditions then imply continuity of the tangential component of the magnetic field as well as of the normal component 
of the magnetic induction when the interface between the superconductor (S) and the magnet (M) is traversed: 

Bt,s = fJ-oHt,M and Bn^s = ^J'^J^oHnM ■ (4) 



A. Infinite flat geometry 



The configuration addressed first is thought to consist of a superconductor extending across the infinite half-space 
—oo < X < Q and a soft magnet extending across the infinite half-space < a; < oo, their interface occupying the 
plane a; = 0, with an externally imposed, homogeneous field Bq pointing in the j/-direction of a cartesian coordinate 
system x, y, z, as shown in Fig. 1. For this usually discussed geometry of the Bean-Livingston barrier™, the presence 
of the magnet does not affect the entry of a straight magnetic vortex parallel to the superconductor/magnet interface. 
Indeed, the vortex self-field here has a tangential component only, which vanishes at the interface owing to the 
vortex image field; the magnetization of the magnet thus is preserved, leaving the vortex field the same as that of a 
vortex near the flat surface of a semi-inflnite superconductor facing vacuum. Nevertheless, the situation could change 
in the more realistic case of fluctuation penetration of a magnetic vortex loop (see Fig. 1), since the magnet then 
might experience additional magnetization, with a corresponding interaction energy contributing to the barrier at the 
superconductor/magnet interface. 

The required solution of equations (1) - (4) for the magnetic induction B can be conveniently decomposed according 
to B = B™ + b, where B™ is the Meissner field induced by the external field in the absence of the magnetic vortex 
loop, 

Bl^{x)^Boexp{x/X), (5) 
We represent the latter field and its source by two- 



and b is the asymptotically vanishing field of the loop itself, 
dimensional Fourier integrals of the kind 



-oo +00 



f(r) 
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dky / dkj^'^^''''^ {x)exp[i (kyU + kzz)] . 



(6) 



-oo —oo 



For a vortex loop situated in the plane z = 0, the Fourier transforms of the cartesian components of the vortex field 
inside the superconductor region are 



i^) = [{biY^^ (0) - P^':f'^y2qX')] exp (qx) 

a 

+ ^ J dx'Q(^Y^ix')exp{-q\x + x'\), 



(7) 



5(fe„.fc.)(^) = (O)exp(gx), 
where q = (fc^ -|- with k — {k^. + A:^)^^^; the boundary values herein read 



[ky ,kz ) 



(fc + iJ,q) A 
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M) (0) 



k {k + ixq) A 



with characteristic integrals of the type 

HY''^ = / dx'Q^ll'^^^ ix') exp i-qx') , 



(8) 



(9) 
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FIG. 1: Schematic view of the semi-infinite superconductor (light shading) and the semi-infinite magnet (dark shading), their 
interface coinciding with the plane a' = of a cartesian coordinate system x, y, z. The solid vertical line indicates the insulating 
layer between the superconductor and the magnet. A straight magnetic vortex parallel to the superconductor/magnet interface 
(arrowed full line) and its image (arrowed dashed line) as well as a magnetic vortex loop (arrowed full semicircle) and its image 
(arrowed dashed semicircle), all situated in the plane z = Q, are shown. The direction of the external field Bo is maxked. 



controlled by the shape of the loop and by the maximum distance between the loop and the interface, a. The vortex 
field inside the magnet region allows the representation b = —fj.iJ,QVip through a scalar potential ip, whose Fourier 
transform is given by 

~(fe k ) pt^v-'--) 

^ ^'^ ^"^^ fc(fe + Mg)A^ '"'^"''^- 

We note that, although the magnetic moment of the vortex loop, with its only component in the direction of the 
external field 

a 

my = J dx'Qf"^^ {x') [1 - exp {-x'/\)] , (11) 



does not depend on fi, the self-energy of the loop, being determined by the field prevailing at the vortex core, could 
well be sensitive to the magnetic environment. Yet, adopting a semicircular shape of the loop with radius a, for which 



g(^fe.,fe.) (3,) = 2i^osm(^kyVa^^^j , (12) 

Ki^) = -2^0 ,f — cos (kyy'a'^ -xA , 
Va^ — a;^ ^ / 

we find, using equations (7)-(10) in the region outside the vortex core whose radial extent may be delineated by the 
coherence length of the superconductor bulk, ^ -C A, the self-energy is represented by the dominant term proportional 
to the length of the loop. 
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apart from small corrections including the factor due to the contribution of the magnet; a result which restates 
that for a semi- infinite superconductor with a flat surface facing vacuumS*^. Intuitively, the decaying tendency of the 
omitted corrections with increasing relative permeability here can be conceived in the following way. The requirement 
of continuity of the normal component of the magnetic induction across the superconductor/magnet interface on the 
one hand, and the definition of the total magnetic flux through its quantization in the superconductor on the other 
hand, ensure the strength of the magnetic induction in the magnet region is typically B ~ $o/A^, whereas the strength 
of the magnetic field in the magnet region is typically H ^ $o//x/iQA^. This yields a contribution to the self-energy 
of the loop proportional to the product of both quantities, which falls of as indicated above. 

The Bean-Livingston barrier arises from a competition between attraction of the vortex loop to the superconduc- 
tor/magnet interface, accounted for by equation (13), and repulsion due to the Lorentz force exerted - by the Meissner 
current - on the loop. In the geometry of Fig. 1, this current flows perpendicular to the plane z = 0, and the work 
done by the external field during growth of the loop from radius a = ^ to radius a ^ A is proportional to the area 
finally covered by the loop, 

^Wf, (Bo) = l^.na^jT (Bo) , (14) 

where 



jT {Bo) = Bo/a^oA (15) 

means the density of the Meissner current at the flat superconductor/magnet interface. From equations (13) and (14), 
the Gibbs free energy of the loop, i.e. the thermodynamic function of relevance here, when ^ <C a <C A, becomes 

^^"^ ~ (iw?) ^^'^ (i) " ' ^^^^ 

identifying the interface barrier against entry of the loop as a function of the strength of the external field. Once, 
when Bq is fixed, the radius of the growing loop has reached its critical size Uc defined by the condition dGfi/da = 0, 
further loop expansion becomes irreversible and vortex entry proceeds. Depending on the quality of the superconduc- 
tor/magnet interface, this may happen at different values of the critical loop radius throughout the range (and even 
beyond) where equation (16) applies. Whilst for an ideal interface, with scale of roughness cr < ^, vortex entry occurs 
at a distance ac = ^, in the case of a real interface, with scale of roughness a < A, vortex entry occurs for Oc = cr. 
Equation (16) in conjunction with equation (15) thus yields for the critical field of first penetration of magnetic flux 
across the flat superconductor/magnet interface 

a form which assumes values between Bd and B^, the lower and, respectively, thermodynamic critical fielcSS,, when 
a varies between A and f . Obviously, the interface barrier against entry of the loop, and therefore the critical field 
of first penetration of magnetic flux as well as the average critical current density deriving from it, are insensitive to 
the magnetic environment, since in the infinite configuration of Fig. 1, the external field remains totally unshielded. 
However, in any finite superconductor/magnet heterostructure, with the range of its magnetic constituent extending 
below the distance to the sources of this field, a significant shielding effect may indeed occur. 



B. Finite curved geometry 



The configuration addressed next is thought to consist of a cylindrical superconducting filament of radius R, 
extended infinitely in the z-direction of a cartesian coordinate system x,y,z and covered with a coaxial magnetic 
sheath of thickness d, the external field Bq again being aligned parallel to the y-direction, as depicted in Fig. 2. 
Whereas the self-energy of an (almost) semicircular vortex loop of radius a <C A, i?, located in the plane and nucleated 
at the filament/magnet interface, duplicates the dominant term given by equation (13), 

with corrections due to the effect of the magnetic sheath even minor than those for the infinite configuration examined 
above^^*^, the Meissner current, and hence the work done by the external field during growth of the loop, as stated 



FIG. 2: Cross-sectional view of the superconducting filament of radius R (light shading) and the coaxial magnetic sheath 
of thickness d (dark shading), their axes coinciding with the z-axis of a cartesian coordinate system x,y,z. The solid ring 
indicates the insulating layer between the superconducting filament and the magnetic sheath. A magnetic vortex loop (arrowed 
semicircle, size not to scale!) situated in the plane z = and nucleated at the filament/magnet interface is shown. The direction 
of the external field Bo and the definition of cylindrical polar coordinates (r, ip, z) are marked. 



by equation (14), does change markedly when considering the finite magnetic environment. We therefore quote the 
Meissner solution of equation (1) expressed in cylindrical polar coordinates (r, Lp, z) adapted to the filament. The 
radial and azimuthal components of the magnetic induction inside the superconductor region are^ 



with 



BoAs [lo (r/A) 
BoAs [lo (r/A) 



4^l 



h [r/X)] sin(^, 
h {r/Xj\ cos If, 



[ii + lf -{^i-lf m/ {R + df IoiR/X) + {n-lf l-Ry{R + d)' hiR/X) 



(19) 



(20) 



where Iq and I2 denote modified Bessel functions of the first kind. Equations (19) in conjunction with Ampere's law 
yield for the density of the Meissner current flowing along the filament. 



jT (r, = 2jT (Bo) Ash (r/X) cos ^; 



(21) 



an expression which adopts its maximum absolute value on the circumference of the filament, Ji™ax(^o)i a-t angles 
ip = and TT indicating the points of most probable nucleation of the vortex loop. Accordingly, the work done by the 
external field during growth of the loop from radius a = ^ to radius a ^ A is 



1 



AW,, (i?o)= 2*o7ra^J-nax(So), 
where, in the practically important limit i? ^ A, 

J^,AB,)^jT{Bo)/a 



(22) 



(23) 
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represents the density of the Meissner current at the curved filament/magnet interface, the parameter 



"^4 



' {R + df 



(24) 



subsuming the combined shielding effect of the superconducting filament and the magnetic environment. Thus, from 
equations (f8) and (23), the Gibbs free energy of the loop, say at the point (p — 0, when ^ ^ a ^ A, becomes 



*0 1^ 2,™ 



Gc. (Bo) = [-^^^ ) ™ [-^J- 2 Jmax (Bo) , (25) 

Evidently, the interface barrier against entry of the loop, taken as a function of the strength of the external field, does 
prove sensitive to the presence of the magnet owing to appreciable shielding of this field in the finite configuration of 
Fig. 2. Again, referring to the condition dGd/da — applied to equation (25) in conjunction with equations (15), 
(17) and (23), we find for the critical field of first penetration of magnetic fiux across the curved filament/magnet 
interface, if i? ^ A, 

Bp = aB°, (26) 

revealing a substantial enhancement by the factor a as compared to the respective field of first penetration of magnetic 
fiux across the flat superconductor/magnet interface for moderately high values of the relative permeability already, 
brought about by the shielding effect. 

The flow of a transport current, of magnitude if, along the filament means that the Meissner current density, given 
by equation (21), is to be supplemented with the corresponding isotropic current densitj*^ 

. , . it Io{r/X) 
'-'^'^ = 2nRXh{R/Xr 

and hence the maximum current density j™ax (^o) entering equation (25) must be replaced by the maximum total 
current density expressed, if i? 3> A, by 

Jtot(Bo)=ra.(So) + ^. (28) 

When Oc is fixed, the condition dGd/da = determines the average critical current density of loss-free transport 
along the filament, jc = it/iTR^, as well. Resorting to equations (23), (25) and (28) in conjunction with equations 
(15) and (17), we get for Bq < Bp, if i? > A, 

Jc {Ba) = , (29) 

where 

B^,^ (Bo) = Bo/a, (30) 

according to equations (19), refiects the maximum strength of the magnetic induction at the curved filament/magnet 
interface. The average critical current density thus is seen to fall off linearly with the strength of the external field, 
at a rate determined by the parameter a, i.e. by the shielding effect of both the superconducting filament and the 
magnetic environment, revealing a strong reduction of the field dependence for moderately high values of the relative 
permeability already, while the zero-field value of the average critical current density is conserved. 

We comment that, if the filament were absent, and hence shielding confined to the magnetic sheath alone, the 
maximum strength of the magnetic induction at the curved inner surface of the sheath would be decreased. By 
formally letting A — s- cx) for fixed R, equations (19) yield 

B'^,^{Bo)^Bo/P, (31) 

with the parameter 



1 

4^1 



{R + df 



(32) 
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duplicating an otherwise derived resulliaS. Since P > a holds for any value of the geometrical and material char- 
acteristics involved, i?°iax (^o) < B^^^^ (Bq) ensues throughout, which confirms that the predicted enhancement of 
the critical field, disclosed by equation (26), like the concomitant attenuation of the external field, revealed by equa- 
tion (30), cannot be exclusively ascribed to the shielding effect of the magnetic sheath, as argued in some previous 
attcmpta^ii4i,42i43^ 

To appraise the relevance of the above results, we take a MgB2 filament with radius R = 5.0 • 10^^ m covered 
by a Fe sheath of thickness d = 2.5 • 10~^ m and relative permeability /i = 50 (RefsiSSi2iS4) , noting the practically 
interesting temperature of 32 K, at which the London penetration depth and the coherence length adopt the respective 
values A = 1.8 • 10~^ m and ^ = 6.5 • 10^^ m (Ref.*"'). If the scale of roughness a, and hence the critical loop 
radius a^, varies between the limits A and ^, the critical field of first penetration of magnetic flux across the curved 
filament/magnet interface, Bp, given by equation (26), is found to range between about 0.16 and 1.02 T (as compared 
to the range between about 0.01 and 0.07 T when the magnetic sheath were absent), and the average critical current 
density, jc (-Bq), from equation (29), turns out to vary between about 6.84 • 10^ and 4.42 • 10** Am-^ at zero external 
field, its rate of change with the field, djc/dBg, amounting to about —4.36 • 10^ Am^^T^^ (as opposed to about 
—6.37 • 10^ Am~^T~^ when the magnetic sheath were absent). These estimates are in accord with the observations 
of recent experiments^^'^'*'^"'. We add that, considering the moderately large ratio of A/^ around 28, the low critical 
temperature of about 40 K and the just minor anisotropy of polycrystalline MgB2 (Refj^^), thermally activated 
penetration of magnetic flux across the barrier, considered befora^2iS, here is insignificant. 

III. SUMMARY 

In conclusion, we have studied the Bean-Livingston barrier at the interface of type-II superconductor/soft-magnet 
heterostructures and demonstrated a characteristic dependence on the geometry of the particular structure and its 
interface as well as on the relative permeability of the involved magnetic constituent. Thus, for the fiat interface 
between a semi-infinite superconductor and a semi-infinite magnet, the external field remains totally unshielded, 
leaving the barrier essentially the same as that at the flat surface of a semi-infinite superconductor facing vacuum. 
However, in any superconductor/magnet heterostructure, where substantial shielding of the external field occurs, the 
modification of the barrier by the presence of the magnet can be significant, as demonstrated for the example of a 
cylindrical superconducting filament covered with a coaxial magnetic sheath. In this finite geometry, with its curved 
superconductor/magnet interface, using typical values of the relative permeability, we predict the critical field of first 
penetration of magnetic flux is strongly enhanced and, concomitantly, the variation of the average critical current 
density of loss-free transport of electric charge with the external fleld is strongly depressed; the zero-field critical 
current density value, however, is retained, since the transport current self-induced magnetic field remains unshielded 
in this geometry. Owing to the expulsion of magnetic fiux out of the filament, the attenuation of the external field, 
and hence the field dependence of the average critical current density, cannot be ascribed to the shielding effect of 
the magnetic environment alone. 
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